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Abstract 

In this paper, we investigate the possibility of explaining both the Rb ex- 
cess and the R c deficit reported by the LEP experiments through Z-Z' mixing 
effects. We have constructed a set of models consistent with a restrictive set 
of principles: unification of the Standard Model (SM) gauge couplings, vector- 
like additional matter, and couplings which are both generation-independent 
and leptophobic. These models are anomaly-free, perturbative up to the GUT 
scale, and contain realistic mass spectra. Out of this class of models, we find 
three explicit realizations which fit the LEP data to a far better extent than 
the unmodified SM or MSSM and satisfy all other phenomenological constraints 
which we have investigated. One realization, the rj- model coming from Eq, is 
particularly attractive, arising naturally from geometrical compactifications of 
heterotic string theory. This conclusion depends crucially on the inclusion of a 
U(l) kinetic mixing term, whose value is correctly predicted by renormalization 
group running in the Eq model given one discrete choice of spectra. 
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1 Introduction & Principles 



During the past six years the four experiments at LEP have provided an abundance 
of data supporting the Standard Model (SM) of particle physics and its SU(3) C x 
SU(2)l x U(l)y gauge group structure. Until recently there has been no significant 
deviation pointing to new sources of physics beyond the SM. However, within the 
last two years there has been growing evidence that a discrepancy exists between 
the predicted and measured widths for the b and c-quark decays of the Z boson. In 
particular, LEP has reported measurements of |TJ: 

R b \ T(Z -» bb/cc) i 0.2219 ±0.0017 ( 
R c }~ T{Z -> hadrons) ~ \ 0.1543 ± 0.0074 Uj 

These values differ from the SM predictions, R b = 0.2152 ± 0.0005 and R c = 0.1714 ± 
0.0001 @ (for m t = (176 ± 13) GeV@ and a s = 0.125 ± 0.010), by 3.9a and -2.3a 
respectively. 

If one is willing to accept the R c discrepancy as statistical, then there are many new 
sources of physics which can serve to resolve the R b measurement by only changing the 
couplings of the third-generation fermions. Such a method is naturally provided by 
low-energy supersymmetry (SUSY) with light charginos and stops or by additional 
fermions mixing with, or additional interactions of, the b and t quarks ||. However, 
if one interprets the R c deficit as another signal of new physics, then the scenarios 
for new physics are more limited |5[]. 

A potential hurdle which one must face with respect to simultaneously explaining 
the R b excess and the R c deficit is that the LEP measurement for the total hadronic 
width of the Z is in good agreement with the SM prediction (Phad = (1744.8±3.0) MeV 
at LEP versus r had = (1743.5 ± 3.1) MeV in the SM), while the sum R b + R c is in 
slight disagreement with the SM prediction. That is, R b + R c = 0.3762 ± 0.0070 as 
measured at LEP (with the error correlations properly included), versus a theoretical 
expectation of 0.3866 ± 0.0005, 1.5<j apart. 

A clue to solving this conundrum may lie in a simple observation. Defining APj 
as the difference between the experimental and the theoretical determinations of Tj, 
one notes that 

3Ar 6 + 2Ar c = (-23.2 ± 24.3) MeV (2) 

so that at the la level, a consistent interpretation of the data is given by assuming a 
flavor- dependent but generation-independent shift in the hadronic Z-couplings. That 
is, 

T u , c = Tl M c +AT c 

r d>s>b = r^ 6 + Ar 6 . (3) 

Such a pattern of shifts has also been suggested in J7|, |[ [|. 
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A second hurdle in explaining the Rb and R c puzzles is that unlike the partial 
hadronic widths of the Z, the well-measured partial leptonic widths are in good 
agreement with the SM predictions: T e = 83.93 ± 0.14 MeV and r inv = 499.9 ± 
2.5 MeV, which are within 0.4<r and —0.4<r respectively of theory. Any source of new 
physics must preserve the successful predictions of the SM for the leptonic widths. 

In this paper we propose to explain the Rb — R c problem by introducing an addi- 
tional U(l)' gauge symmetry. If this new U(l)' is broken near the electroweak scale, 
there can be significant mixing between the usual Z and the new Z'. The physical 
Z-boson as produced at LEP will then have its couplings to fermions altered by an 
amount proportional to the Z — Z' mixing angle times the Z' coupling to those same 
fermions. 

Analyses have recently appeared in the literature || ||| that seek to fit the LEP 
data by introducing such an additional U(l) f . Both of these works make a phenomeno- 
logical fit to the data introducing some number of new parameters, such as arbitrary 
U(l)' charge ratios, Z — Z' mixing angle, and Z 1 mass. These analyses do indicate 
that this class of scenarios has the potential to solve the Rb — R c discrepancy, and 
are therefore interesting. However, they share some fundamental problems associated 
with the lack of an underlying, consistent framework. For example, the extra U(l)' 
is not anomaly free (this is true both for the [£7(1)'] 3 , and most seriously, the mixed 
SM-[/(l)' anomalies). Further, since the authors of 0, |9[ also seek to explain the 
CDF dijet excess, they are forced to take a high value of the Z' mass. For such Z' 
masses, the f/(l)'-couplings have to be so large that the U (1)' gauge coupling becomes 
non-perturbative at most a decade above the Z' mass scale; implicit in this is that 
the Z' width in these models equals or even exceeds the Z' mass. 

Here we will take a different approach. We set forth a few basic principles which 
we believe any attractive Z'-model should obey. Within this framework we will find 
that there exist only limited classes of U(l)' models which are phenomenologically 
viable and theoretically consistent. Each class has a well-defined prediction for the 
U(l)' charges of the SM fermions, reducing much of the arbitrariness in the couplings. 
We will not attempt to explain the CDF dijet anomaly. 

The principles that we demand are: 

• The low energy spectrum must be consistent with the unification of the standard 
model gauge couplings that occurs in the minimal supersymmetric standard 
model (MSSM). This will lead us to consider models which are extensions of the 
MSSM, with any non-MSSM matter added in particular combinations which can 
be thought of as filling complete multiplets of SU(5). We allow the possibility 
of unification within a string framework, and do not require the presence of a 
field theoretic GUT. 

• All non-MSSM matter must fall into vector-like representations under the SM 
gauge groups. Such a requirement is consistent with the absence of experimental 
evidence for new fermions with masses below the top quark mass. Further, 
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note that additional chiral matter is disfavored by the electroweak precision 
measurements, since, in contrast to vector-like matter it can give very large 
contributions to the 5, T, and U parameters. 

• The U(l)' charges of the SM leptons must be (to a good approximation) zero. 
This requirement of leptophobia is motivated by the phenomenology. This alone 
will eliminate the U(l) factors associated with most traditional GUT groups, 
since GUT's tend to place leptons and quarks into common multiplets. 

• Consistent with Eq. (|3|), we require that the U(l)' couplings be generation- 
independent. This requirement is essential if tree-level hadronic flavor changing 
neutral current processes mediated by the U(l)' gauge boson are to be naturally 
suppressed. This also has the advantage of simplicity and economy. 

To be precise, the principle of unification that we will impose requires that the 
meeting of the SM couplings at 2 x 10 16 GeV is not a coincidence. For simplicity we 
will not explicitly consider in this article the various string models where the scale 
of unification is increased to the (weak-coupling prediction of the) string unification 
scale Mgfr loop ~ 5 x 10 17 GeV, such as those discussed in IfLOfl , although it will be 
clear that the consequences for our discussion of such a modification are slight. (Note 
that one interesting possibility that could maintain unification at 2 x 10 16 GeV is the 
strongly coupled string scenario recently proposed by Witten JO] . ) 

If one takes the unification of gauge couplings to imply the existence of a simple 
GUT gauge group, then the natural candidates with extra £/(l)'s and three chiral 
families are SO (10) and Eq. However the single additional U(l) within 50(10) is not 
leptophobic. In Eq all linear combinations of the two additional C/(l)'s orthogonal 
to hypercharge couple to leptons. Nonetheless, we will show that by including an 
effect usually overlooked in the literature (L r (l)-mixing in the kinetic terms through 
renormalization group flow [|T^, |I3"| ) there exists a unique U(l)' in the Eq group which 
is compatible with the data. The Eq subgroup in question is usually known in the 
literature as the ?7-model and interestingly is the unique model which results from Eq 
Wilson-line breaking directly to a rank-5 subgroup in a string context |14j]. We will 
discuss this case in some detail in Section f|. 

Finally, although we assume the MSSM for the purposes of gauge coupling unifi- 
cation, we do not use MSSM loop contributions to the Zbb vertex in order to explain 
any part of the anomaly. In particular we do not assume light charginos or top 
squarks which are the necessary ingredients for such a scenario H. 



2 Z -Z' Mixing 

We begin with a brief general discussion of Z — Z' mixing in the context of an 
SU (2) L x U(1)y x U (1)' model. A more detailed discussion can be found, for example, 
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in Refs. [15[ |16 |. The neutral current Lagrangian of the Z and Z' is given by 



^nc = \ E Al" (—(Vi + a a 5 )Z, + g\v[ + a'rf)Z'\ ^ (4) 

where 

Vi = T 3i — 2QiS w , cii = — T 3 j (5) 

are the SM vector and axial couplings of the Z, and v', a' are the (unknown) vector 
and axial couplings of the Z'. Here g' is the coupling constant of the new U(l)' and 
Sw = sin 2 9 W . 

After electroweak and U(l)' breaking, the Z and Z' gauge bosons mix to form the 
mass eigenstates Z\^ where we will identify the Z\ with the gauge boson produced 
at LEP: 

Z\ = cos £ Z + sin £ Z' 

Z 2 = — sin£ Z + cos£ Z'. (6) 

Since such mixing must necessarily be small in order to explain the general agreement 
between LEP results and the SM, we will throughout this paper use the approximation 
Z\ ~ Z + £Z'. We will also assume that the mass of the Z 2 is large enough so that its 
effects at LEP, either via direct production or loop effects can be ignored. Therefore all 
new physics effects must appear through the mixing angle £. The relevant Lagrangian 
probed at LEP will then be 

£zx = 7T~ E ^i-ffr + W 5 )Zi^i (7) 

ZCw j 

where, for small £, 

Vi ~ Vi + 1v[ 

a>i ~ Oi+la't, (8) 
and we have defined the auxiliary quantity 

i={g'c w /g 2 )fi. (9) 
Because the Z\ is no longer purely the electroweak Z, the p-parameter 

p - 1 = aV2 g f (n u (o) - n 33 (o)) (10) 

receives a tree-level correction. (Here Iljj(O) are the SU(2)l vacuum polarization 
amplitudes at zero momentum transfer.) If we define the corrections to p by 

p EE 1 + Ap SM + Ap, (11) 
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where Ap$M is due to loop corrections already present in the SM (such as the top), 
then the mixing with the Z' contributes to Ap. Since we will later be interested in 
taking into account the effects of further shifts in p due to the rest of the MSSM 
spectrum, we decompose Ap = Ap M + Ap extra , where Ap M is the part due to mixing 
with the Z' . The value of Ap is the quantity that our fits to the LEP data will directly 
constrain. Writing the Z — Z' mass matrix as 

- ( 5 % ) ■ w 

then for M§, 3> m z , Am 2 , one finds that the shift in p due to mixing, ApM, is given 
by 

Ap „^(^W(MR (13) 



where 

Am 2 



Ml, 



(14) 



There is also a corresponding shift in s 2 ^: 



z 

2 



■ 2 — a 2 I S W C W A n 

"w ~ s w 



In terms of the above parameters, one can then calculate the Z\ partial width to 
fermions: 

T(Z 1 ^Jf) = ^^pN c (v} + a}). (16) 

A further relation may be obtained by examining the specific form of the terms 
that come into Eq. (|T2|). If we assume that the fields fa which receive vev's occur 
only in doublets or singlets of SU(2) L , then 

C W i zc w 
M z , = 2g' 2 Y,(Q^) 2 , (17) 

i 

C W i 

where Q\ is the U(l)' charge of 0, and v\ is the sum of the vev's of the SU(2)i 
doublets. Then we may write Ap^j as a simple function of £: 

a pm^-(Jt-) {^)l=- A 4T,^){Q'A). (is) 

\9 cw J \ m z J v z . 
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What is noteworthy about this relationship is that it is connects the two quantities 
(ApM and £) which are experimentally constrained at LEP (up to Ap ex tra, which we 
can bound), in a way that is independent of the unknown gauge coupling g' and the 
Z' mass. Note that Am 2 and £ in Eq. flHf ) have opposite signs, so that Apu is always 
positive. 



2.1 U(l) a x U(l) b Mixing and RGE's 

The discussion so far has echoed the conventional wisdom on the subject of Z — Z' 



mixing. However, it was realized many years ago JT^| that in a theory with two 
£7(1) factors, there can appear in the Lagrangian a term consistent with all gauge 
symmetries which mixes the two £7(1) 's. In the basis in which the interaction terms 
have the canonical form, the pure gauge part of the Lagrangian for an arbitrary 
U(l) a x U(l)b theory can be written 

r — ^ T?^ v W, > ^ T?V V TP, , S * n ^ rpiiv TT, , 

+ Am 2 A ia) ^ b) + -m 2 a A {a) ^ a) + -mlA^Afc (19) 

If both [7(l)'s arise from the breaking of some simple group G — *> U(l) a x U(l)b, then 
sin% = at tree level. However, if the matter of the effective low-energy supersym- 
metric theory is such that 



E (QiQi) * o, (20) 



i=chiral fields 



then non-zero \ will be generated at one-loop. This is naturally the case when split 
multiplets of the original non-Abelian gauge symmetry, such as the Higgs doublets 
in a grand unified theory, are present in the effective theory. Since we are interested 
in a large separation of scales, M GUT and M z , we will need to resum the large log- 
arithms that appear |13|, [l7j using the renormalization group equations (RGE's) for 
the evolution of the gauge couplings including the off-diagonal terms. 

Once a non-zero x ( or Am 2 ) has been induced, one needs to transform to the mass 
eigenstate basis. To do so, one must perform a (non-unitary) transformation on the 
original gauge fields, A^ and Aq), to arrive at the mass eigenstates, Z 12 : 

A(a) = (cos£ — tan%sin£) Z\ — (sin£ + tanxcos£) Zi 

A(b) = (sin£Z 1 + cos£Z 2 ) /cosx, (21) 

where 

— 2 cos y (Am 2 — m 2 sin y) . . 

tan2£ = — = ° - AJ . (22) 

mi — m l a cos 2x + 2 Am 2 sin \ 

This transformation results in a shift in the effective charge to which one of the 
original Z7(l)'s couples. (One U(l) can always be chosen to have unshifted charges.) 
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This can be seen by taking the £ = limit of the above transformation. The resulting 
interaction Lagragian is then of the form 1 12] : 



Ant = *l> Y {QaQaZlf, + (g b Qb + g a bQa)Z 2 ^) ij) (23) 

where the redefined gauge couplings are related to the original couplings, g°, by g a = 
0a > 9b — g b / cos X an d g a b = — 0a tan X- The ratio 5 = g a b/gb is a phenomenologically 
useful parameter, representing the shift in the i^-fermion coupling due to kinetic 
mixing. 

The renormalization group equations for the coupling-constant flow of a U(l) a x 
U(l)b theory, including the off-diagonal mixing, are most usefully formulated in the 
basis of Eq. (E3f). In this basis the equations for the couplings g a , g b and g a b are: 



^ 1 ,3p 

dt " 16tt 2 ^ aa ' 

dgb 1 / o 2 - \ 

~dt = m^ 9b * 9b + 9ab + 2gbgabBab I ' ( 24 ) 
(it 16vr 2 



g'tgabBbb + flr^Soa + 2g z a g ab B aa + 2g' z a g b B ab + 2g b g 2 ab B, 



at) 



where S^- = ti (QiQj) with the trace taken over all the chiral superfields in the 
effective theory, and there is no sum over (a, b) in Eq. (0). From these equations 
we immediately see that even if g ab = to begin with, a non-zero value of the off- 
diagonal coupling is generated if the inner-product ti(QiQj) between the two charges 
is non-zero. The advantage of this basis for the RGE's is that the low-energy value of 
the parameter 5 is given directly by the ratio g a b/g b evaluated at the low scale. (This 
is not the case for the more symmetrical form of the RGE's given in Ref. |l3|j .) 

For the case at hand, we will choose the couplings of the usual to be canonical, 
shifting the charge of the Z' . Since it is the B^ component of which mixes through 
the kinetic terms, the couplings of the Z' to matter fields can be expressed in terms 
of an effective U(l)' charge Q e fr = Q' + YS, where Y is hypercharge. We can translate 
from Eq. (|23f) using g a = gilcyj and g b = g' so that g a b = — #2 tan#vf tanx and 
^ = gab/gb- The vector and axial couplings that come into Eq. (|8|) are given by 

v 1 = Q cS W-Q cS (iJ c ) 

a' = -Q eS (ip)-Q eS (ip c )- (25) 

Note that both if) and ip c are left-handed chiral fields: Q e s(ip c ) = —Q c fi(ipR)- 

In most of the models we will consider, we will work directly with Q e s; in such 
models, whether or not Q e ff can be expressed as some Q' + Y5 for non-zero S will 
not have an effect on the analysis. However, when considering the 77-model coming 
from Eq, the difference between Q e ff and Q v will have important consequences on the 
observable physics. We reserve further comment on the U(l) mixing in the Eq model 
until Section El 
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Kinetic mixing of C/(l)'s will also shift the p-parameter. In the previous subsection 
we had assumed that we could write the electroweak Z in terms of the mass eigenstates 
as Z = cost; Z\ — sin£Z 2 . However, in the presence of a non-zero x (or 5), this is 
changed to (see Eq. (pT|) , replacing tanx with — s^ytanx): 



Z = (cos£ + sin£svy tanx) Z\ — (sin£ — cos^s^y tanx) Z 2 

Z 1 = (sin£ Z\ + cos£ Z 2 )/ cosx (26) 

A = 7 — cy^ tan x(sin£Zi + cos £Z 2 ) 

where 7 is the physical photon. Eq. (|22|) for £ becomes 

, of -2 cos x( Am 2 + m|siy sin x) , , 

tan zcj — — ~ ~ - — J - — - , v J 

Mf, — cos^ x + rn-zSw s ^ n X ~ 2 Am^s^/ sin x 

while the Zi mass is given to lowest order in mfz/Mf, by: 

9 2 I m l / Am 2 \ 2 

m^^m^l-^^ + ^sinxj >. (28) 

The coefficient of the Z\ term in Eq. (^) is essentially a wave-function renormal- 
ization for the Z\ and contributes to ApM by absorbing part of the explicit mass shift 
which came from mass matrix mixing |16| . The net effect is a negative contribution 
to Apj\/ which subtracts from the positive definite contribution coming from mass 
mixing. In terms of 5, 

M 2 

Ap M ~ - 2k Z 5 ( 29 ) 

where k = g'cwsw / 92- The important point to note is that, in the presence of kinetic 
mixing, ApM can be smaller than had there been no such mixing; in fact, Ap^j can 
be negative. 

The kinetic mixing also shifts s 2 ^ beyond what was already included in Eq. ([To]): 

s 2 w = s 2 w k=s=o - & ( + ks) . (30) 

For 5 = this reduces to Eq. flT5|). Finally, there is a new contribution, Sm, to 
the so-called S-parameter (see, e.g., Ref. [|16|]) due to kinetic mixing which can be 
negative: 

aS M ^ -4c 2 w k£5 (31) 

to leading order in ml/Mf,. 
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2.2 New Contributions to Oblique Parameters 

As noted in the previous Sections, in the absence of U(l) kinetic mixing (i.e., 5 = 0), 
Z-Z' mixing gives a positive contribution to the p-parameter, denoted by ApM, and 
no contribution to the ^-parameter. Since our numerical fits are sensitive to the total 
Ap and S, it is important to see if there are corrections from sources other than the 
Z-Z' mixing. (Both Ap and S are defined to be zero in the SM for some reference 
top quark and Higgs boson masses which we take to be 175 GeV and 125 GeV re- 
spectively.) The spectrum of the effective theory in all models that we will consider 
includes a Higgs sector with two doublets, vector-like states in complete U SU(5) mul- 
tiplets" , and the superpartners of all particles, all of which can in principle contribute 
to the oblique parameters. The sizes of these contributions depends on the details 
of the mass spectrum. As we shall see, the scale of the U(l)' breaking turns out 
to be relatively low in all models (typically Mz> ~ 200 — 250 GeV). Therefore the 
contributions of the additional matter cannot be ignored in general. Let us therefore 
estimate the typical allowed ranges for Ap ex t ra and S ex t V3u (S = Sm -\- S ex t ra j, given 
some reasonable choices for the spectrum, in particular that depending upon MSSM 
superpartners, Higgs sector, and additional vector-like matter. 

The superpartner contributions to Ap extra and S exkTa in the MSSM have been 
studied in Refs. and [[Lj| respectively. In Ref. (T^j it has been shown that such 



contributions to S cxtra _ are generally very small; therefore we will ignore MSSM su- 
perpartner contributions to S^xtra i n everything that follows. Likewise it is shown in 
Ref. |Tj| that the corrections to Ap extra from the MSSM sparticle spectrum are small 
(and positive) with the exception of the stop-sbottom correction which can be sizable 
depending on the nature of the supersymmetric spectrum. 

Although the Higgs-boson contribution to Ap extra in a general two-doublet model 
can be large and negative (as large as —0.01), in supersymmetric models there are 
restrictions on the Higgs sector parameters, resulting in an absolute lower bound of 
Apextra > —0.0015 from the MSSM Higgs sector. However, in the class of models 
which we will consider in Section [|, this number becomes —0.002 since the Higgs 
sector in these models is not identical to that of the MSSM. This is because the 
pH u Hd term of the MSSM will be replaced by XH u HdS, where S is a SM-singlet 
field carrying U(l)' charge. There is also a new contribution to the Higgs potential 
from the U(l)' D-term. We have analyzed the Higgs spectrum of these models, which 
resemble the MSSM with a singlet (the NMSSM). In the limit where the singlet vev 
is large compared to the doublet vev's, but keeping the mass of the pseudoscalar 
fixed, we have numerically examined the most negative Ap cxtra obtainable from the 
Higgs sector and found it to be —0.002. Of course, this could be partially offset by 
some positive contribution from other sectors, such as the stop-sbottom sector. In 



the model analysis of Section |3.1| we will therefore consider two cases, one in which 
Ap extra = and another in which we take Ap extra to have the not unreasonable value 
-0.001. 

As far as the contributions from additional vector-like matter are concerned, we 
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will always consider the simple isospin-symmetric case (i.e., the masses of the T 3 = ±~ 
states equal) where there are no vector-like contributions to Ap extra . In this limit, 
Sextra need not be zero. For the various models we will consider, S extr & receives 
potentially large contributions from the multiplicity of lepton/higgsino doublets which 
arise. There are two natural cases. One, in which the vector-like contributions to the 
doublet masses dominate over the chiral contributions, gives S extTau ~ 0. Alternatively, 
because the weak scale and the U(l)' scale are quite close, the chiral masses can be 
of order the vector-like masses; we have estimated, using the results of Ref . |2(J , the 
contribution to S^tra in this case to be +0.14 per pair of such doublets. 



3 Leptophobic U(l) Models 

Any model which hopes to extend the SM in a minimal fashion must give masses 
to the SM fermions through the usual Higgs mechanism. Within a supersymmetric 
model, such couplings appear in the superpotential, W. Letting Wq be the minimal 
superpotential consistent with the SM, we writeS 

W = h u QH u u c + h d QH d d c + h e LH d e c . (32) 

The new U(l)' must also preserve this superpotential. Demanding that the U(l)' 
couplings of the leptons be zero allows us to write the charges of the remaining fields 
as: 

Q'(Q) = x Q'(H U ) = -x-y 

Q'(u c ) = y Q'{H d ) = (33) 

Q\d c ) = -x 

We next require that the resulting gauge theory have no anomalies. In the case 
of the SM particle content alone, this implies C 3 = C 2 = C\ = C = 0, where, 



[SU(3)} 2 x 17(1)' 
[SU(2)} 2 x 17(1)' 
[U(l)v] 2 x U(l)' 

[U(l)f X U(l)y 



3x + 3y = C 3 (34) 
8x - y =C 2 (35) 



-x + ly = C x (36) 
(x + y)(7x-5y) = C . (37) 



At this time we do not concern ourselves with the [U (l)'] 3 , or U(l)' [gravity] 2 anomalies 
since these can be saturated with any number of SM gauge singlets. The only solution 
which cancels all anomalies in Eqs. (|34|)-(|3~?|) is the trivial solution x = y = 0. 

Going beyond the MSSM, we wish to add matter in such a way that the unification 
of gauge couplings that occurs in the MSSM is not upset. To do so we must arrange 
that the additional matter changes the MSSM one-loop beta-function coefficients in 



1 With the extended matter content that we will introduce later in the paper, it is also possible to 
consider more complicated non-minimal choices for these Yukawa couplings, where the Higgs that 
couples to e c and d c are distinct. We will not analyze these possibilities in detail here. 
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such a way that A6 2 = A6 3 = §A&i. This constraint can be most easily understood 
as requiring the addition of complete 577(5) multiplets to the spectrum (though U(l)' 
need not commute with this fictitious SU(5)). 

Our principles outlined in Section [l] constrain us further in how we add SU(5) 
multiplets to the model. Implicit in the requirement of unification is that the gauge 
couplings remain perturbative up to the unification scale. This implies that we can 
only add (a limited number of) 5's, 10's, and their conjugate representations. By 
requiring that all new matter be vector-like under the SM gauge groups, we restrict 
ourselves further to adding the multiplets in pairs. In combination, these two princi- 
ples limit us to adding (A) up to four (5 + 5) pairs, or (B) one (10 + 10) pair, or (C) 
one pair each of (5 + 5) and (10 + 10). 

Consider Model A with a single pair of (5 + 5). Because we require neither that the 
U(l)' commutes with the ersatz 577(5), nor that the charge assignments be vectorial 
with respect to the U(l)', we write general U(l)' charges for the new states as: 

5 = (3,1) [-1/3, a x ] + (1,2) [1/2, a 2 ] 

5 = (5,1) [1/3, + (1,2) [-1/2, a 2 ] ^ 

where each state is listed by its (577(3) c , SU(2)l) [U(1)y,U(1)'] representation/charge. 
The anomaly coefficients are changed to: 

Cq — > Cq — a\ + a 2 + a 1 ~ a 2 C 2 — > C 2 + a 2 + °2 /oqn 

C\ — > G\ + |(ai + ai) + |(a 2 + a 2 ) C 3 — > C 3 + a x + a 1 . 

Solving for the condition C 3 = C 2 = C\ = C = yields 

y = 2x, (40) 

with the additional relations a\ = — 2(a 2 +9x)/3, a 2 = — a 2 — 6x, andai = (2a 2 — 9x)/3. 
Note that all charges are rationally related, and, further, that for a purely axial choice 
of U(l)' charges {a\ = a,\ etc.), the only solution is the trivial one x = y = a, = 0. 
The result Eq. ( [40"D does not depend on the number of (5 + 5) pairs. Thus for this 
entire class of models, we know the couplings of all the quarks to the Z' through 
Eq. (|33"D, up to one overall normalization. 

The same exercise can be undertaken for Model B. Now we add the states in the 
(10 + 10) with charge assignments 

10 = (3,2)[l/6,a 3 ] + (3,1) [-2/3, a 4 ] + (M)[l,a 5 ] Ul] 
10 = (3,2)[-l/6,a 3 ] + (3,1) [2/3, a 4 ] + (1, 1) [-1, a 5 ]. 1 1 

In the general case the phenomenologically important ratio y/x is undetermined by 
the anomaly conditions. However, if we make the very natural simplifying assumption 
that the Z7(l)' charges in Eq. (f4T|) are purely axial (a 3 = a 3 , etc.), then the [f/(l)'] 2 x 
U(l)y anomaly equation (p7|) is unmodified and there are only two solutions for the 
charge ratio: 

7x 

y = -x, or y = —. (42) 
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The associated charges of the extra states are {a 3 , a 4 , a 5 } = {— 3x/2, 3x, — 3x/2} and 
{— lla;/10, — 7x/5,x/10} respectively. In the following we will refer to these models 
as "B(-l)" and "B(7/5)". In the "B(-l)" model the charges are identical to baryon 
number, with the Higgs doublet H u carrying zero charge. At this stage it is important 
to recognize that both these models have the potential problem that the extra states 
do not include (1, 2)±i/2 representations which can be used to give a naturally small 
off-diagonal mixing term Am 2 in the Mf z , mass matrix Eq. (|T2"|). In the B(-l) 
model, there is no tree-level Z — Z' mixing. Even at the one-loop level, no such 
mixing arises in the simplest version of this model where the (10 + 10) states receive 
masses from SM singlets only. In the B(7/5) model, on the other hand, there is 
tree-level Z — Z' mixing, which however tends to be too large. As we will see, this 
model requires additional (negative) contributions to the p-parameter to relax the 
constraint Eq. (|18|) . 

Model C has, in the general case, ten new U(l)' charges corresponding to the 
ten new states in Eqs. (|38|) and (fill) , and again even with the constraints imposed 
by anomaly cancellation the ratio y/x is not determined. However there are two 
particularly attractive and natural subclasses of these models. In the first subclass the 
U(l)' charges of the extra states are chosen to be purely axial. This leads to the charge 
ratios y/x = -1 or 7/5 as in Eq. © (Models "C(-l)" and "C(7/5)" respectively). 
Note that since all C-type models contain an extra pair of Higgs doublets, they are 
naturally able to accommodate a suitably small Z — Z' mixing. The second attractive 
subclass of Model C is defined by setting the U(l)' charges of the anti-generation 
(5 + 10) to zero (a x = a 2 = a 3 = a 4 = a 5 = 0). In this case the ratio y/x is 
continuously adjustable as is the charge, a 3 , of the additional (3, 2)i/ 6 state. Among 
this continuous family, the choice 

y = x (43) 

is especially simple and attractive (Model "C(l)"). 

In all cases we still need to impose the [[/(l)'] 3 and U(l)' [gravity] 2 anomaly can- 
cellation conditions. It is important to consider an efficient way of achieving this 
because we will soon see that there is a strong constraint arising from the require- 
ment of perturbativity of the [/(l)'-coupling all the way up to the GUT scale, and 
the U(l)' beta-function gets a significant contribution from these SM-singlet states 
(collectively E's). One must also add sufficient vector-like states charged under U(l)' 
to give all the additional matter (including states both in the 10 + 10 and 5 + 5's, 
and the E's) masses. The derivation of the minimal set (in the sense of reducing 
their contribution to the beta-function) of states and charges that satisfies these 
conditions is a difficult problem in general. As our interest is only in the value of 
the minimal U(l)' beta-function coefficient b (including the contributions from the 
SM-non-singlets states) we just quote the results for 6 min for the various models in 
Table 1, and where we have employed an ansatz for the spectrum of anomaly can- 
celling states!. (Our ansatz is to choose a set of U(l) '-charged states, S, which cancel 

2 We doubt that it is possible for some of the SM singlets to be very light, which would have 
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Model 


A B(-l) B(7/5) C(-l) C(7/5) C(l) 


^min 


1363 280 174 129 154 191 



Table 1: Minimal beta- function coefficients (in the normalization x = 1) for the mod- 
els defined in the text, together with additional SM-singlet matter to cancel [C/(l)'] 3 and 
gravitational anomalies, and give mass to all non-MSSM states. The version of Model A 
considered has a single 5 + 5. 



the extra anomalies and simultaneously contribute minimally to the U(l)' /3-f unction. 
We then include a minimal set U(l)' vector-like states which give mass to the S's.) 

Strictly speaking our "unification principle" does not absolutely require the per- 
tubativity of U(l)' up to the GUT scale - it is only the SM gauge couplings that 
we require to successfully unify while still perturbative. For instance, it is possible 
that our extra U(l)' gauge symmetry is enhanced into a non-Abelian gauge symmetry 
well before the GUT scale, in which case the following is (possibly much) too severe 
a restriction. Nevertheless it is interesting to see the bounds on the mass of the Z' 
that follow from such a requirement. 

The restriction is derived as follows: Using the Eqs. (^) and ( |T3D for the fitted 
quantities £ and ApM, we find that for the x — 1 normalization choice, 

However requiring that the Landau pole does not occur until a scale A gives (at one 
loop) the restriction 

^ * Ti^W (45) 

where b is the beta-function coefficient. Putting these two equations together leads 
to a restriction on the Z' to Z mass ratio in terms of the "measured" quantities £ 
and Ap^j, and the coefficient b (for which we have a lower bound given the minimal 
spectrum of U(l)' charged particles necessary for anomaly cancellation, etc.): 

fMzV <142 ApM 1 ( } 

\M Z J ~ (0 2 6bg(A/M z )' 1 ! 

For the most restrictive case of A = 2 x 10 16 GeV, this gives 

AfeV £ 4 3 



reduced significantly the 0- function coefficients 6 m in- Constraints on this possibility come predomi- 
nantly from supernova cooling and to a lesser extent big-bang nucleosynthesis (BBN). If these SM 
singlets are massless, they will be produced copiously inside supernovae through their Z' interactions. 
Once produced, they will free stream out of the supernova leading to rapid cooling. Consistency 
with SN1987A observation requires that the Z 1 mass must be greater than about 1 TeV or that the 
singlet states must be heavier than about 30 MeV. 
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o 


7T^\ 

4(1) 

s± 


7rT\ 

/\ 


rn — 

B 


77T\ 

Uq 




0.98 


-1.02 


-0.55 


0.50 


Re 


-0.04 


-0.83 


-0.78 


0.71 


O'had 


0.006 


0.12 


0.32 


-0.29 


Rb 


0.007 


0.16 


-2.8 


-0.71 


R c 


-0.004 


0.33 


5.4 


1.4 


M w /M z 
A b 

FB 


0.38 


-1.0 











-56 


-2.1 





A c 
FB 





-59 


2.4 


-5.4 


At 

rt-FB 





-115 









Table 2: Coefficients Aq and Bq and observables O used in the fit to the electroweak data, 
as defined in Eqs. (pi)- (|4£ 



3.1 Experimental Constraints 

Having defined each class of models, we know that each will, by definition, be lepto- 
phobic. However it remains to be seen if they can describe the physics as observed at 
LEP any better than the SM. Note that as far as the agreement with the LEP data is 
concerned, the only important feature of a model is the value of the ratio y/x. (In all 
models except the ?pmodel of Section |] we will choose to normalize the U(l)' gauge 
coupling g' such that the quark doublet charge x = 1.) 

To study this question, we have performed a \ 2 fit of each model to the LEP 
data, broadly following the procedure of Refs. | 
observables as inputs: T z , Re = Thad/r^, cr had , R b 
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R r 



We take 9 independent LEP 
M w jM z Ab 



-c "*w/mz, A b FB , A C FB , and A FB . 
Theoretically, the shift in each observable O can be expressed as a function of Ap, £, 
x, and y: 

— = A a Ap + (Bg>x + B$y) I (48) 

However, it is only in the simple case of no kinetic-mixing that expressions for Aq 
and Bq follow directly from those given in Refs. || 0. This is because they take 
Eq. ( |T5| ) as the relation between s 2 ^ and Ap M ; that is, the expressions of Refs. || pT5| 
assume that 5 = 0. For 5^0, Eq. ( j30|) holds instead. We then re-express 



A a Ap 



AffAp + A^As 2 



w ■ 



A 2 — 2 _ 2 



w\i=&=o 



(49) 



where A 



(i) 



implicit dependence through As^. The coefficients Aq 
the discussion of Ref. IT 



; numerical values for the A}n and B%' are given in Table 2. 



iq includes only the explicit dependence of the observable O on Ap, not the 

ri,,, .„., ,(r;,-.;,„. + . \ 1 '■) are eas jiy generalized from 

Note that As^ is not a new parameter to be fit, since it is simply a function of ApM, 
£ and 6 through Eqs. ( p9"D and (P0|). Clearly for 5 = the procedure here reduces to 
that of Refs. g [TJ|. 

Unlike Ref. M, we have opted against using the data from SLC. As is well known, 
the SLC data is approximately 2a from the corresponding data at LEP. This could 
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Model 


Ap 






x 2 


a s (M z ) 


SM 


5 x 1(T 5 







22.8 


0.125 


2 


9.1 x 10~ 4 


-4.6 x 10" 


-3 


10.9 


0.125 


-1 


-5.6 x 10~ 4 


-4.1 x 10" 


-3 


14.8 


0.110 


7/5 


3.5 x 1(T 4 


-7.6 x 10" 


-3 


5.4 


0.125 


+ 1 


-2.6 x 10" 4 


-8.9 x 10" 


-3 


4.0 


0.123 



Table 3: Results of fit to LEP data in the Standard Model (at a s = 0.125, the best fit for 
the LEP data alone) and models with charge ratios y/x = 2, —1, 7/5, +1. In all cases the 
X 2 are for 7 dof, and m t = 175 GeV and m-Higgs = 120 GeV are assumed. The best fit 
value of a s in the range 0.110 to 0.125 is quoted in each case. 



be a systematic effect at LEP, SLC (or both), or a sign of new physics. Here we will 
take this discrepancy not to be a sign of new physics. Therefore, as the effects we are 
studying (Rb and R c ) are in the LEP data, we choose, in this paper, to exclude the 
SLC data from our fits. 

In our fits for the models of this Section, we have taken S ex t ISL = and allowed for 



Ap cxtra to be either zero or —0.001 consistent with our discussion in Section [2.2| . The 



negative value of Ap extra in particular leads to a relaxation of the mass limits on the 
Z'. 

In Table 3 we have shown the x 2 f° r each of the possible charge ratios y/x = 
2, —1, 7/5, and +1 in addition to the SM; the SM is defined by setting £ = in the 
fit. For each model, we have given the values of Ap and £ at the minimum x 2 , as well 
as the value of a s in the range 0.110 < a s < 0.125 which produces the best fit to the 
data. For two of the models listed, the best fit value of Ap is negative; however, the 
fit depends only weakly on Ap so that positive values of Ap are allowed at relatively 
low x 2 as shown in Figure 1. 

For the two most attractive models, C(7/5) and C(l), we have included plots in 
Figures 2 and 3 of iso-x 2 contours in the (£, Apjv/) plane. The solid ellipses represent 
contours of x 2 = 14.1 and 18.5, values which correspond to goodness-of-fits of 95% 
and 99% respectively for 7 dof, assuming Ap extra = 0. In both cases, the contours 
impinge significantly into the physical Apm > region. The dashed ellipses represent 
the case for which Ap extra = —0.001 as discussed earlier in the text; for this case the 
allowed values of Ap M are larger. 

Figures 2 and 3 also show contours of constant M z > calculated assuming the 
perturbativity constraints of Eq. (|47|) and using the values of b m i n tabulated in Ta- 
ble 1. For the C(7/5) model, the 95% (99%) C.L. bound on M z > is 180 (350) GeV for 
Ap cxtra = and 250 (500) GeV for Ap oxtra = -0.001. Similarly, for the C(l) model 
the 95% (99%) C.L. bound on M z , is 150 (300) GeV for Ap cxtra = and 220 (450) GeV 
for Ap extra = —0.001. The B(7/5) model has mass limits only slightly stronger than 
those of the C(7/5) model: 170 (320) GeV for Ap extra = 0. For the remaining models 
in Table 1, the corresponding Z' mass limits are much stronger (with the exception of 
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the ^-model of Section f|, which falls into the broad class of model A but has smaller 
value for the /^-function coefficient b). 

One might expect that Z'-models of the type considered here would be strongly 
constrained by either UA2 or CDF/DO. However, the strongest Z' mass bounds 
in the literature depend on observation of the leptonic decays of the Z', which are 
highly suppressed in these leptophobic models. The dijet decays of the Z', which 
dominate its width, are hard to detect above background except for limited ranges 
of Z' masses and couplings. In particular, CDF can only exclude Z' — > jj for Mz> 
roughly between 400 and 460 GeV [^l]|, and then only for SM strength (or stronger) 



couplings. UA2 has a similar bound of Mz> > 260 GeV |22| , but here again one 
requires SM strength couplings. Note that because of the small couplings that result 
from our perturbativity constraint, we tend to find that the production cross-section 
for the Z' at a hadron collider is suppressed by at least 40% compared to the SM Z 
cross-section. We therefore find that UA2 does not provide a strong constraint on 
the Z' mass in these models. 

All of the theoretical mass bounds that we have derived depend strongly on the 
value of the U(l)' gauge coupling, and thus on the size of b min and especially on 
the assumption of perturbativity of the U(l)' gauge coupling all the way up to the 
GUT scale. If the U(l)' interaction is enhanced to a non-Abelian group at some 
intermediate scale, then the Z' mass bounds are much weaker; we are investigating 
this possibility. By either decreasing 6 min or decreasing A (the scale up to which 
we require perturbativity), g'(M z ) will increase. As g' increases the Z' mass bound 
increases but the Z' production cross-section at a hadron collider, relative to a Z of 
the same mass, also increases. At some mass, however, the kinematic suppression of 
the Z' production wins and the experimental bound goes away. We will not consider 
the details of these competing effects here. 

Taking all the phenomenology together, including the possibility of naturally small 
Z-Z' mixing, we view the C(l), C(7/5), and the //-model of the next Section as 
promising Z' explanations of the R^, R c anomalies. 



4 The 77-model 

As we noted in Section [j], E® is a natural, and for our purposes, minimal, choice for a 
simple GUT group containing extra £7(1) 's. In addition E§ appears as an underlying 
feature in many geometric compactifications of the E§ x Eg heterotic string. In either 
case, the list of possible subgroups into which the Eq can break is small and well- 
defined. 

Since E 6 is rank-6, its Cartan subalgebra contains two U(l) generators besides 
those of the SM gauge groups. At scales just above the electroweak scale, the addi- 
tional gauge symmetry could appear either as a commuting U(l)' factor (as we have 
been assuming up to this point) or as a unification of the SM groups into some non- 
Abelian group (e.g., SU(4) C x SU(2)l x SU(2)r). The latter choice cannot describe 
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Figure 1: 99% C.L. contours for the four basic classes of models labeled by their Q/u c 
charge ratio in the (£, Ap) plane. The cross represents the SM. 



the physics at LEP since it cannot be leptophobic. Returning to the former, we can 
write the new U(l)' as a combination of the two extra t/(l)'s in E e , usually denoted 
as U{l) x and U{1)^\ 

Q'(a) = cosaQ x + sinaQ^. (50) 

In Table 4 the charges Q x and Q^p are given for each of the states of the MSSM using 
the standard embedding into the 27. 

No linear combination of U(l) x and U{1)^ is completely leptophobic. The best 
one can do is to find models for which the axial coupling of the charged leptons 
is zero. Since the vectorial contributions for charged leptons appear proportional 
to 1 — ~ 0.07, the Z' coupling to charged leptons could be highly suppressed 
with respect to the hadronic couplings. However, such models would necessarily have 
couplings to the neutrinos of order the hadronic couplings. If, after Z-Z 1 mixing the 
net effect were an increase in r inv at LEP, the model could be quickly ruled out. On 
the other hand, if r; nv were to decrease, one could imagine that some new source of 
invisible Z-decays {e.g., neutralinos) could offset the difference. We consider such a 
scenario to be fine tuned and do not consider it here. 



However, as was discussed in Section |2.1| , in an arbitrary U{l) a x U(l)b model, 
there is one more free parameter, a mixing parameter g ab for the two groups. In 
the case of the breaking of some unified gauge group, Ggut, at some high scale into 
Ggut -> SU(3) C x SU(2) L x U(1) Y x U(l)', the value of g ab will be zero at the high 
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Figure 2: x 2 contours for the C(7/5) Model in the (£,A/9a/) plane. The solid ellipses 
represent the 95% and 99% C.L. bounds on the fit. The dashed ellipses represent the 
corresponding bounds if Ap cxtra = —0.001. The three solid lines are contours of Mz' arising 
from the theoretical constraint of perturbativity of the U (1)' coupling up to the GUT scale, 
and are labeled in GeV. 
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Table 4: U(l) charges of the states of a 27 of Eq. 
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Figure 3: \ 2 contours for the C(l) Model in the (£, A/j/y) plane. See caption of Figure 2 
for explanation. 

scale. Nonetheless, through its RGE's, Eq. (|4]), g a b will be driven to non-zero values 
for generic particle content. The effective coupling to the Z' is then Q eff = Q'(a) + SY 
where 5 = g ab /g'. 

From the low-energy point of view, 5 is a completely free parameter which must 
be fit to the data just as we did £ or Ap. Therefore, we have repeated the x 2 analysis 
of the previous Section; however the charges of the SM fermions are now completely 
determined in terms of a instead of x and y. Figure 4 is a x 2 plot in the plane of (a, 5) 
showing the fits to the LEP data at 95% and 99% C.L. At each point in the plane, the 
X 2 value is minimized with respect to the remaining two free parameters, Ap and £. 
Along the bottom of the plot are indicated the values of a consistent with the Xi ip, 
and i] models (a = 0, ir/2, — tan -1 ^5/3 ~ —0.91 respectively) commonly discussed 
in the literature. All previous discussions of these models (with the exception of 
Ref. |23|) have tacitly taken 5 = 0. 

What is remarkable about the fit is that it picks a very particular model out, for 
a limited range of 5. To fall within the 95% C.L. region (x 2 < 14.1), a model must 
have a = -0.89 ± 0.06 and 5 = 0.35 ± 0.08. Recall that the SM has af = 22.8 in 
the same parameterization. Only one model lies within the region of allowed a: the 
so-called ?7-model. The charges of the MSSM states under U(l) v are given in Table 4. 

That the best fit in the (a, 5) plane lies at Q' ~ Q v and 8 ~ 1/3 is not surprising. 
The effective charge Q e s = Q v + Y/3 is completely leptophobic; in fact it is the only 
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a 

Figure 4: x 2 contours for general Eq models. The two contours represent confidence levels 
of 95% and 99%. Three canonical Eq models are labeled at the bottom. The two points 
highlight the remodel with 5 = 1/3 (x) and 5 = 0.29 (A). 

combination of the three Abelian generators in Eq which is leptophobid. Note that 
the Qrj charges of the lepton doublet L and the lepton singlet e c are proportional to 
their hypercharges. Thus, U(l) v is uniquely picked out as capable of describing the 
new physics at LEP. In Figure 4 we have shown the 5=1/3 rj- model with a cross. 

If U(l)' is indeed U(l) v , there are a number of direct consequences both for theory 
and phenomenology. First, U(l) v does not fit into any GUT group smaller than Eq. 
Thus, if the unification of the gauge couplings at a scale near 10 16 GeV is not an 
accident, it indicates either a true field-theoretic Eq GUT (and no SU(5) or 50(10) 
unification) or string-type unification in which SU(3) x SU (2) l X U(l)y x U (l) v unifies 
directly at the scale Mmssm = 2 x 10 16 GeV. Second, cancellation of the anomalies in 
Eqs. (l34f)-(l37D requires the existence of three complete 27's of Eq. Besides the usual 
states of the MSSM, one can expect three pairs of D and D c quarks which are SU(2)l 
singlets with Y = =Fl/3, two additional pairs of SU{2) L doublets with Y = ±1/2, 
three right-handed neutrinos, u^, plus SM singlets (at least one of which will receive 
a vev to break U(l) v and will be eaten by the Z'). 

We can now write the mass matrix of the Z-Z' system. Denning tan (3 = (H u ) / (Hd) 

3 After submission of this paper, we were kindly informed by F. del Aguila that the possibility 
of a leptophobic U(l) in Eq had been observed in Ref. |2?j]; however, it was not realized that the 
required value of 8 was naturally generated through radiative effects in a model with realistic matter 
content. 
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and g v to be E^-normalized, the off-diagonal element in the mass matrix is given as 
in Eq. flTJ): 

Am 2 = ^l^iTsi^dQv + SY)^) 




2c w U 9n92V2zSin2(3 (51) 

where the last equality holds for the case where the only SU(2)l doublets with non- 
zero vev's are H u and H^. For the completely leptophobic //-model (i.e., 5 = 1/3), £ 
and Ap M are then simply 



t 9v c w /5 . 2n(M 2 z \ 5. ( 1 



Unfortunately, such a relationship between ApM and £ does not provide a very good 
fit to the data except near the unphysical value of tan/3 ~ 0.6; the best fit consistent 
with Eq. (|2|) and tan/3 > 1 has \ 2 of 22.0, not much better than the SM \ 2 = 22.8. 
There is a second related problem: since Am 2 ~ Mf and we expect (in the absence of 
tuning) for the Z' mass to be only somewhat heavier, we should expect large mixing 
angles £ to result. This is generic problem of U(l)' models where the U(l)' is expected 
to be radiatively broken close to the weak scale [J24| . 

The solution to both problems involves the introduction of additional SU(2)l 
doublets, charged under U(l)', which receive vev's near the weak scale. In our case 
these will play several roles: arranging the /3-functions of the model to unify at the 
GUT scale, allowing for small £ by cancelling the H u contribution to Am 2 , likewise 
decoupling Ap M from £, and driving 6 > 0. 

Consider, for example, extending the minimal mmodel to include the pair of dou- 
blets which fit into the [78, 16 + 16, 5 + 5] of [E 6 , SO (10), SU(5)}, with the doublet 
in the 5 getting a vev, vi, near the weak scale. Then in the leptophobic 77-model, 
Am 2 oc (t>| sin 2 /? — vf). If a near cancellation can be arranged between the two 
terms in Am 2 , then small mixing will result and simultaneously ApM <C £ as needed 
phenomenologically. Since Mf oc (v\ + vf) and we need vz and vi of the same order, 
the Higgs vevs, v u and Vd, which give masses to the fermions will be proportionally 
smaller. In the case Vd <^ v u ~ Vi, the large top-bottom mass ratio is natural and 
the top Yukawa is of the same size as one would expect in the MSSM with tan/? = 1. 
This is actually still below the top Yukawa infrared pseudo-fixed point, which now 
takes a larger value (/if xed ~ 1.25) because of the slow running of a s in this model. 

Imposing on the superpotential of the minimal mmodel a discrete Z2 symmetry 
(a simple extension of the usual i?-parity) one finds: 

W v = Qu c H u + Qd c H d + Le c H d + SH u H d + SDD C + Lv c E u (53) 

Under the i?-parity, all the states of the 27 are odd except H u , H d and S. This su- 
perpotential forbids dimension-4 proton decay; dimension-5 operators are also known 
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to be unobservably small in the //-model [25]. There appears in the superpotential a 
Yukawa mass term for the right-handed neutrino fields, Lv c H u . To be consistent with 
current neutrino mass bounds, this coupling must be small or zero or the v c must 
have large Majorana mass terms through some singlets. By flipping the .R-parity as- 
signment of the v c one can forbid the term altogether, but at the price of introducing 
into the superpotential the term v c Dd c . Such a term would lead to D-d c mixing were 
v c to receive a non-zero vev. 

One can also expect radiative symmetry breaking much as in the MSSM. If the 
SDD C coupling is 0(1), the soft mass term for the S'-field, m|, will be driven negative 
through its RGE's, triggering [7(1 ^-breaking through (S) ^ at a scale just above 
the electroweak scale. (The electroweak symmetry will similarly be broken by m 2 Hu 
running negative due to the large top Yukawa coupling.) Since the singlet S has no 
electroweak interactions unlike H u , it is conceivable that the mass-squared of the S 
fields turns negative at a larger momentum scale compared to H u . The non-zero (S) 
will in turn produce a [iH u Hd and a fi'DD c term. For SH u Hd and SDD C couplings 
of 0(1), one expects ~ Mz>- In particular, it is natural for the D and D c states 
to be heavier than the Z. Finally, we note that there is no mechanism within the re- 
model for v c to receive a vev radiatively which does not violate some other constraint 
(such as neutrino mass bounds) f25| . Thus D-d c mixing will not occur. 

The //-model with only three 27's of E 6 does not satisfy all of our initial principles 
because it does not have gauge coupling unification. As mentioned above, unifica- 
tion can be arranged by introducing one pair of SU(2) l doublets with hypercharges 



5/3 Qy = ±|- From a string point of view, these may be viewed as coming from 
a 27 + 27 or a 78, the rest of whose states received masses at the string scale 
This, along with anomaly cancellation considerations, requires the doublets to have 
equal and opposite Q v . If these doublets also have non-zero effective charges Q v + 5 Y, 
their vev's may contribute to the Z-Z' mixing matrix as outlined above. A problem 
may potentially arise in trying to generate vev's for these doublets radiatively; one 
possibility is to allow couplings of the type H u H' d through singlets. 

(This model has, beyond the spectrum of the MSSM, three each of (3, 1) and (3, 1) 
and six of (1, 2). This is exactly the content of three (5 + 5)'s of SU(5). Note that in 
terms of the charge ratio y/x, the purely leptophobic (5 = 1/3) //-model is equivalent 
to Model A of Section [3|. However, the presence of kinetic mixing (5 ^ 0) induces 
contributions to the oblique electroweak parameters not present in Model A. Also 
unlike the purely leptophobic models of that Section the value of 5 in the ?7-model is 
generically not 1/3, but is instead determined through the RGE's and thus through 
the low-energy spectrum. Further, its /3-function is substantially smaller than that 
of Model A with a single (5 + 5), since for the //-model the anomaly cancellation is 
generation by generation, providing a more economical set of charges.) 

There are two variants of the //-model for which the value of 5 at the electroweak 
scale is of particular interest: (i) The "minimal" //-model that possesses three gener- 
ations of 27's and one additional vector-like pair of Higgs doublets that arises from 
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Table 5: Beta-function coefficients for the Minimal and Maximal 77-models, GUT normal- 
ized. 



the 78 of E 6 . These doublets have charges y5/3Qy = —1/2 an d 2^/l5Q v = 6 under 
the GUT-normalized U(1)y x U(1) v symmetries; (ii) The "maximal" 77-model with 
in addition to the states of the minimal 77-model a further effective 5 + 5 of SU(5) is 
added (so that unification is preserved) , but which is composed of a second vector-like 
pair of the doublets in the 78 together with the color triplets D + D coming from 
the 27 + 27. The maximal model has the largest field content consistent with pertur- 
bative unification of the gauge couplings at 2 x 10 16 GeV. The values of the charge 
inner products Bij for these two models are given in Table 5. The field content of 
both these models is consistent with small Am 2 in the Z-Z' mass matrix. 

Running the SM couplings up to the unification point and then numerically run- 
ning the RGE's of Eq. (|4]) for gy, g v and gy v down to the electroweak scale, we find 
predictions for 5 in the two models: 

5 min = 0.11, 5 max = 0.29. (54) 

Both of these are calculated with a s (Mz) = 0.120. Larger values of a s (Mz) lead to a 
slight increase in the values of 5 compared to Eq. (|54|). The threshold corrections to 
5 coming from mass splitting of the light states are typically of order 0.01. It is quite 
remarkable that the totally leptophobic value of 5 = 1/3 is very nearly predicted by 
the renormalization group running of the "maximal" ?7-model. From the one-loop 
RGE's, the value of the U(l) v gauge coupling at the electroweak scale is g v = 0.40. 

Given these values of 5 we can now investigate how well the ?7-model variants 
can fit the LEP data. As discussed in Section [2.2| we will consider both the case 
of S'extra = and Sextra = 0-14 P er P an ~ °f higgsino/lepton-like doublets. We will 
take Ap extra = 0. The minimal model is clearly disfavored by the data, having a 
X 2 no better than the SM for both values of S^xtra- Likewise the maximal model 
with Sextra = is disfavored. The phenomenologically favored maximal model has 5 
doublet pairs giving S^xtra — 0.7 and a minimum \ 2 — 13.9 at a Z' mass of 215 GeV; 
this is within the 95% C.L. bounds shown in Figure 4, where the model is indicated 
by a triangle. At the minimum, S = Sm + S'extra = —0.1. Note that the goodness of 
the fit does not depend strongly on the exact value of S^tra in the range 0.5 to 1.5; 
in particular the resulting S only varies within the range —0.1 to 0.1. 

Given g v and the bounds on Ap M and £ we are in a position to calculate the bounds 
on the Z' mass, using Eq. (|i~3"D. For the //-model with 5 = 0.29, we find that in order 
to fall within the 95% (99%) C.L. limits for our fit, then M z , < 240 (420) GeV, under 
the assumption of no additional contributions to Ap. (New positive contributions to 
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Figure 5: x 2 contours for the 77-model with 5 = 0.29 in the (£, Apm) plane. See caption of 
Figure 2 for explanation. Additional positive contributions to Ap reduce the best fit value 
of the Z' mass. 



Ap, which are natural in these models, push the best fit Z' mass to lower values.) 
These fits are shown in Figure 5. UA2 has performed a Z' search in the dijet channels, 
excluding a Z' with 100% branching fraction to hadrons and SM strength interactions 
up to masses of 260 GeV However, given the value of g v = 0.4 and the U^)^ 
charges of the quarks, one can show that the production cross-section for this Z' is 
approximately 1/4 that of the Z, too small to be excluded at UA2. 

What is remarkable about this analysis is that the ?7-model, which has been ex- 
tensively studied in the literature and for which strong bounds on its mixing with 
the Z and its mass have been published, has been resuscitated by the inclusion of 
the additional U(l) kinetic mixing effect. This is even more so, since the value of 5 
is correctly predicted in specific models in which only one discrete choice of matter 
content has been made! 



5 Conclusions 

In this paper, we have investigated the possibility of explaining the Rb excess - R c 
deficit reported by the LEP experiments through Z-Z' mixing effects. We have con- 
structed a set of models consistent with a restrictive set of principles: unification of 
the SM gauge couplings, vector-like additional matter, and couplings which are both 
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generation-independent and leptophobic. These models are anomaly- free, perturba- 
tive up to the GUT scale, and contain realistic mass spectra. Out of this class of 
models, we find three explicit realizations (the rj, C(7/5), and 0(1) models) which fit 
the LEP data to a far better extent than the unmodified SM or MSSM and satisfy 
all other phenomenological constraints which we have investigated. The //-model is 
particularly attractive, coming naturally from geometrical compactifications of het- 
erotic string theory. This is especially so since the value of the mixing parameter, S, 
is correctly predicted given only one discrete choice of matter content. 

In general, these models predict extra matter below 1 TeV and Z' gauge bosons 
below about 500 GeV, though the Z' of these models will be difficult to detect exper- 
imentally. 

Note Added 

After this work was completed two further interesting works concerning the exper- 
imental consequences of leptophobic £/(l)'s appeared PB|[]2l|. These papers noted 
that there can exist important low-energy constraints on leptophobic models arising 
from atomic parity violation (APV) and deep-inelastic neutrino scattering experi- 
ments. In particular, Ref. |28] argued that the aesthetically appealing models that 
we have constructed in this paper are strongly disfavored by the APV data. While 
this is usually true in the heavy Z' mass approximation that we have been employing 
up to now, this conclusion does not hold in the very interesting case of a light Z' 
(Mf, > as we will now outline. 

The APV experiments result in constraints on the so-called weak nuclear charge 
Qw of various elements such as Cesium and Thallium with high atomic and neutron 
numbers Z and N. The charge Qw is itself defined in terms of the product of 
the axial electron coupling with the up and down type quark vector coupling via 
Q w = -2 {C lu (2Z + N) + C ld {Z + 2N)} where 

G 

v 2 i=u4 

In the case where the Mz> ~ %, both Z x and Z 2 exchange contribute to the coeffi- 
cients Cu. In the approximation where the mixing is small (< 1, but no expansion 
is made in the mass ratio r = (m|/Mf/), the expression for the Cus is 

Cu = -{vi+vg(l-r)}, 

where the v 's and £ are defined in Eqs. (|J) and (§) respectively. It is therefore clear 
that the constraint from the APV data becomes vacuous as r — > 1. Specifically, we 
find that the APV data do not significantly increase the total x 2 f° r Z' masses below 
about 150 GeV. 
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One may similarly consider the effect of a leptophobic Z' on the neutrino scattering 



experiments. We find that the parameters e l L and e l R defined in Ref. ||30|| , are altered 
by an amount 

respectively. Thus the weaker constraints from the neutrino scattering data also 
dissappear for light to moderate Z' masses. 

We will address the full fit including these constraints (as well as the SLC and 
other data), more fully in a forthcoming paper Ref. pi|, where we will also discuss 



the models with variant Higgs structure mentioned in the Section [| footnote. 
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